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Abstract 


In Nekrasov and Shatashvili pointed out that the M = 2 instanton partition function in 
a special limit of the Q-deformation parameters is characterized by certain thermodynamic 
Bethe ansatz (TBA) like equations. In this work we present an explicit derivation of this fact 
as well as generalizations to quiver gauge theories. To do so we combine various techniques 
like the iterated Mayer expansion, the method of expansion by regions, and the path integral 
tricks for non-perturbative summation. The TBA equations derived entirely within gauge 
theory have been proposed to encode the spectrum of a large class of quantum integrable 
systems. We hope that the derivation presented in this paper elucidates further this com- 
pletely new point of view on the origin, as well as on the structure, of TBA equations in 
integrable models. 
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1 Introduction 


A number of remarkable connections have been observed between gauge theories and inte- 
grable systems. They appear to be useful to increase our understanding of both subjects. On 
the one hand, using powerful integrability techniques one may hope to solve certain gauge 
theories non-perturbatively. On the other hand, gauge theory can help to formulate and 
solve integrable system. A spectacular example is given by planar M = 4 super Yang-Mills 
theory, for which in the last ten years or so tremendous progress has been achieved in solving 
the theory based on the underlying integrability and on the AdS/CFT correspondence, see 
[I] for a recent review. There exists another amazing connection between gauge theories and 
integrable models. The gauge theories in this case are M = 2 supersymmetric gauge theo- 
ries. They do not need to be planar but the connection with integrable models is restricted 
to a special class of supersymmetric observables. In this paper we focus on an important 
object in this class, the so-called instanton partition function and its relation with quantum 
integrable systems. 


In the groundbreaking work of Seiberg and Witten [2] [3], the exact solution for the low 
energy effective action of certain NV = 2 gauge theories was proposed based on holomorphicity 
properties and electromagnetic duality. The low energy dynamics are encoded in a single 
object, called prepotential F(@). It is a holomorphic function on the Coulomb moduli space, 
with coordinates @, and can be reconstructed from the so-called Seiberg-Witten (SW) curve 
and SW differential. Shortly after, it was realized that this description provides a direct 
connection between M = 2 gauge theories and classical algebraically integrable systems 
[4] 56], see e.g. [7] for a pedagogical introduction. These are essentially a complex analogue 
of integrable systems in the sense of Liouville. 


The challenging program of obtaining the Seiberg-Witten prepotential by a direct gauge 
theory calculation, developed on [8] 9| [10], was finalized in [LI]. This problem was solved us- 
ing powerful localization techniques. Interestingly, this calculation produced a two-parameter, 
called €; and €, deformation of the prepotential. The SW prepotential can be obtained as 


F(a;q) = lim «ee. log Z(đ@, 61, €2;¢), (1.1) 
€1,€2—0 
where we add the explicit dependence on the coupling constant q, but suppress dependence 
on further parameters such as masses. The partition function Z receives tree level, one loop 
and instanton contributions. The latter part is usually referred to as Nekrasov instanton 
partition function. 


The parameters €1, €2 correspond to a Lorentz rotation Q, thus the name ()-deformation, 
that encodes certain twisted boundary conditions for the four dimensional gauge theory. It 
was first introduced in [8] [9] in order to regularize the volume of the instanton moduli space. 
The Q-deformation can be understood in a simple way by considering the five dimensional 
lift of the M = 2 theory, further compactified on a circle [12]. In this set-up it is inter- 
preted as twisting of R* by a Lorentz rotation while going around the circle. Introducing 
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a two-parameter generalization of the prepotential triggered a huge progress. An impor- 
tant example is the connection with topological strings. Upon taking €& = —€2 = gs, the 
gauge theory partition function reproduces the topological string partition function with gs 
as genus parameter [14]. The question of what is the topological string theory analog 
of the €, Æ €) case leads to the definition of so-called “refined” topological strings [16]. 
Another spectacular example of progress driven by the calculation of the Nekrasov partition 
function is given by the Alday-Gaiotto-Tachikawa (AGT) correspondence [7]. 


More recently, Nekrasov and Shatashvili proposed that upon taking the limit €&2 > 0 
and interpreting € as Plank constant, one obtains a correspondence between supersymmetric 
vacua of a given gauge theory and eigenstates of the corresponding quantum integrable 
model. The relation between the SW prepotential and classical integrable systems is thus 
quantized. This is usually called Nekrasov-Shatashvili (NS) limit and will be the main focus 
of this paper. The central role in this correspondence is played by the so-called twisted 
superpotential W(đ, €1; q) defined as 


WG, €1;q) = lim ez log Z(ä, €1, €2; q) - (1.2) 
E=} 


The Nekrasov-Shatashvili proposal is that, once this function is known, the eigenstates of the 
quantum integrable system are classified by solutions of the following quantization condition 


exp (=w) =), ere ae (1.3) 


ay 


where r is the rank of the gauge group. These equations identify the twisted superpotential 
W with the so-called Yang-Yang (YY) function of the quantum system. The proposed 
correspondence provides an efficient general mechanism to define and solve quantum integral 
models. Remarkably, it can be argued that the instanton part of the prepotential Winst, 
defined via (1.2), can be characterized as the solution of certain non-linear integral equation 
of the Thermodynamic Bethe Ansatz (TBA) type [20]. The main goal of this paper is to 
develope some of the ideas presented in to give a more explicit derivation, as well as 
some generalizations, of such TBA equations. 


The proposal above originates as some sort of extension of the so-called Bethe/gauge 
correspondence 22]. The latter is based on the observation that the vacuum equations 
of two dimensional M = 4 gauge theories, broken to M = 2 by twisted masses, coincides 
with Bethe equations for integrable models. The two dimensional twisted superpotential 
is equal to the YY function that encodes the Bethe equations and Coulomb parameters 
correspond to the Bethe roots. The generators of chiral ring of the gauge theory are 
mapped to Hamiltonians of the integrable model, while their expectation values mapped to 
the corresponding eigenvalues. In this way one obtaines a large class of integrable models 
whose spectrum is characterized by traditional, possibly nested, Bethe equations. Many in- 
tegrable models do not belong to this class. The simplest example is given by the quantum 
Toda chain, see e.g. [24], whose classical limit is connected to four dimensional pure SU(N) 
SYM. From the insight of the Bethe/gauge correspondence it is then natural to propose 
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the two dimensional twisted superpotential which corresponds to such integrable models. It 
is the effective low energy action for the four dimensional M = 2 gauge theory subject to 
an Q-background that preserves two dimensional M = 2 super-Poincaré symmetry, namely 
€2 = 0. This observation provides a stong motivation for the above correspondence. An 
interpretation of this correspondence was given using brane constructions in [25]. A fur- 
ther essential step in understanding the nature of the relation between quantum integrable 
systems and gauge theories has been presented in [26]. 


The proposal of Nekrasov and Shatashvili has inspired many other studies. Let us men- 
tion a few. In it was shown that, similarly to the prepotential F, the twisted 
superpotential W(e,) can be obtained by calculating period integrals of a suitably deformed 
SW differential. This analysis is also inspired by the AGT correspondence, by which the NS 
limit corresponds to the semiclassical limit of Liouville CFT BOJ. 


In our work the Coulomb parameters will be assumed to be in generic positions. Extra 
considerations are needed if they take special values. For example in the conformal SU(N) 
SYM with Ny = 2N, if the Coulomb parameters a are set to be equal to mfuna — 2 €1, 
where n € Z, one can quantize the corresponding integrable system to obtain a lenght 
N spin-chain with infinite dimensional heighest weigth representations of slp at each site 
[33]. Its spectrum is described in terms of traditional Bethe Ansatz equations. Such 
developments triggered the discovery of a number of new dualities between various integrable 
models [34] [35] B6, [37] [38]. The NS proposal has also inspired various studies in (refined 
topological) string theories [39] [40] [41] [42] [43] [44] [45] [46] where the the general Q-background 


plays a crucial role. 


Despite the importance of this correspondence, the precise mechanism by which the 
instanton part of the twisted superpotential defined in turns out to be characterized as 
the solution of TBA equations is still to be elucidated. In this paper, we will fill this gap. In 
order to fully prove the NS’s proposal, at least in some example, one should be able to show 
that the same TBA equation characterizes the spectrum of the corresponding integrable 
model. In the case of pure SU(N) SYM, which corresponds to the periodic Toda chain with 
N sites, this was achieved in [24]. This interesting problem will be studied elsewhere [47]. 
In the following we briefly outline the main ingredients used in our derivation of the TBA 
equations for (1.2), as well as the structure of the paper. 


As pointed out in [18], it is convenient to start with the contour integral form of the 
instanton partition function. In this representation the instanton partition function can be 
interpreted as the partition function of a non-ideal gas of particles. The particular structure 
of the two-particle interaction potential makes the study of the e2 — 0 limit rather subtle. 
More precisely, this potential is the sum of a short-range (of order €2) strongly attractive 
piece and a long-range interaction part. In Sections P.I] and 2.2] we consider the simplified 
situation in which either the long- or short-range part is set to zero. In order to study 
these simplified partition functions in the €2 > 0 limit, we combine a number of techniques 
like Mayer expansion (a standard method in statistical mechanics) and the method 
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of expansion by regions [49] (a powerful method to compute Feynman integrals in small 
parameter expansions). For the case with only long-range interactions, the €2 — 0 limit 
turns the logarithm of the partition function into a sum over certain tree graphs. On the 
other hand, the free energy corresponding to only short-range interactions gives rise to the 
dilogarithm function Lig, which can be shown either by direct residue calculation of relevant 
integrals or via Mayer expansion together with the method of expansion by regions. 


In order to study the full partition function we find it convenient to use an iterated 
version of Mayer expansion, see [50]. This expansion effectively creates a new partition 
function whose “new particles” are clusters of the original particles. The interaction within 
each cluster is governed by the short-range interaction, the one between different clusters 
by the long-range part. The iterated Mayer expansion thus produces an expression for the 
twisted superpotential Winst as a sum over tree graphs, with vertices given by clusters. This 
expansion is carried over in some details in Section The expression can be compared 
to high order in the instanton number with the expression coming from the solution to the 
TBA equation as discussed in Section [8.1] providing direct non-trivial check of the equality. 


There is an elegant way to prove that this equality holds to all orders in the instanton 
counting parameter q. It is based on rewriting the grand canonical partition function of the 
non-ideal gas in terms of a (0+1)-dimensional path integral. The analysis needs some special 
care as the potential has an unusual feature of depending in a singular way on h, which is 
identified with € in this case. A slight modification of the argument in [5I], together with 
the calculation of the contribution from the short range interactions corresponding to the 
dilogarithm, shows that the instanton partition function in the e2 — 0 limit is obtained by 
the saddle point evaluation of the path integral. The saddle point equations are nothing but 
the TBA equations. This is explained in details in Section 


In Section [4] we present a generalization of the TBA equations corresponding to quiver 
gauge theories. More precisely, we consider quivers characterized by a Dynkin diagram of 
ADE, or ADE type. The twisted superpotential for such theories is shown to satisfy a 
set of coupled TBA equations with one equation for each node of the quiver and couplings 
corresponding to edges in the quiver. The derivation is a simple extension of the one for 
the single gauge group case. This is so as the short range interaction, responsible for the 
clustering of particle, is non-vanishing only for particles corresponding to the same gauge 
group factor in the quiver. 


In order not to overload the main text, in Appendices we include some review material 
together with a few technical points concerning the derivation. A review of the contour 
integral form of the instanton partition function is given in Appendix[A] Some useful formulas 
are collected in Appendix [B] A discussion of the method of expansion by regions is given in 
Appendix [C] In Appendix [D] we present an alternative derivation of a tree graphs expansion 
of the instanton partition function. 


The full partition function is a product of three terms Z = Ztree Z1-loop Zinst- In this 
paper we will be only concerned with the study of the instanton part Zj,,. For this reason 
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from now on it will be simply denoted by Z. 


1.1 NS’s correspondence 


To complete the introduction, we present the integral representation of the instanton parti- 
tion function for pure SU(N) N = 2 super-Yang-Mills and the corresponding TBA equation. 
The derivation of the TBA starting from the gauge theory expression of the instanton par- 
tition function is the main goal of the paper. 


Instanton partition function 


The instanton partition function can be written in a contour integral representation as 


a= S) E (6) T] Dn), (1.4) 


J>I 


where $77 = ¢; — 67, € = € + €2, and 


D(a) = > (1.5) 


1 
The origin of this expression is reviewed in Appendix [A] The parameters entering this inte- 
grals, namely €12 and a, are taken to be real with a small positive imaginary part 70. The 
domain of integration above should be understood either as a real slice integration or equiv- 
alently, upon closing the integration in the upper-half plane, as a multiple contour integral. 
In Appendix[A.4] we review how the residue evaluation of (L4) reproduces the representation 
of the instanton partition function as sum over N-tuples of Young diagrams. 


We emphasize that the precise form of Q(x) and D(a) does not affect the derivation 

presented in this paper. This is the main reason why the generalization to quiver gauge 
2 

theories is rather straightforward. On the other hand the presence of the factor =o in 


2 
D(x), which have a particularly singular limit for e2 small, will play a crucial role and will 
be responsible for the appearance of the dilogarithm function in the TBA. 


TBA form 


The claim of is that the twisted superpotential, defined as 


W = Limit. 40 (€2 log Z), (1.7) 
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can be written as the critical value of the following integral functional: 


1 fdr dy dx ph 
vod =} | HH ece- now + | E poele) + Li Qe") ] a8) 
which is 
AR ~ol@)) 4 Li ~o(z) 
W = Crito V A] = | x5 |- 522) log (1 — qQ(x)e®™) + Liz (GQ(a)Je ) | } (1.9) 
where v(x) is the solution of the following TBA-like equation 
dy —(y) 
(x) = | zg G(x —y) log (1 —aQ(yje"™) (1.10) 
Q(x) is defined in (1.6), and the propagator G(x) is related to D(x) as 
Ee a (z = +) (1.11) 
€9 dx T= Ey 


2 Mayer-Cluster expansion 


As mentioned in [18], the contour integral form of the instanton partition function (4) can 
be interpreted (for each k) as the partition function of a one dimensional non-ideal gas of 
particles ¢),...,¢,, subject to an external potential and a pair-wise interaction potential 
respectively given by 


U(x) = — log (2m) , Vi™(z) = —log (D(z)) . (2.1) 
1€2 

Upon summing over the number of particles k, the instanton partition function takes the 
form of a grand canonical partition function. The free energy of this gas can be studied by 
Mayer expansion techniques (see for example for a nice introduction), as pointed 
out in [18]. In this section, we will perform such kind of expansion in full details. The 
limit of €2 —> 0 appears to be rather subtle. In order to perform this limit we need to face 
the problem of studying the leading behavior of some multiple integral where a parameter 
is small. It turns out that this can be conveniently studied by employing the method of 
expansion by regions [49] discussed in Appendix [C] We introduce this method to provide a 
unified framework to study certain integrals, but we stress that all the result of this section 
are also obtained without employing such technique. 


In order to analyze the behavior of the partition function (1.4) in the limit in which € is 
small, it is convenient to split the function D(x), see (1.5), into two parts: 
Po a = g =e 
D(x) = ——— D(z), Dlr) = =—>. (2.2) 


Dae Dae 
Eo = ey gpi = e 
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The reason of such decomposition is as follows. The factor -e corresponds to a pair-wise 
interaction which is strong and attractive at distances of order €2 and rapidly decreases at 
large distances. The remaining factor D(x) corresponds to a pair-wise interaction which is 
different from zero only at distances of order e. Thus corresponds to splitting the 
potential into short- and long-range parts. The natural candidate to study such kind of 
potentials is the so-called iterated Mayer expansion [50]. As we will see in Section 2.3] this 
effectively creates a new grand canonical partition function whose “new particles” correspond 
to clusters of the original particles. We will start by considering some simplified situation. 


2.1 Only long range interactions 


2 
Let us first consider a simplified version of (1.4) without the factord i, i.e 
EJ 2 


Zig = 2 i Jg Wer) I] Pl). (2.3) 


2 clei<k 
The basic idea of Mayer expansion, see e.g. [52], is to introduce the function fry as 


Dig) = 14+ fry, (2.4) 


and expand the interaction products as 


TL (+f) =1+ fnt D fufes +... (2.5) 


1<I<J<k I<J LTT ZJ" 


Each monomial in the right hand side of this equation can be visualized as a graph (not nec- 
essarily connected) on the set [k] = {1,2,..., k}. More precisely, each particle in {1, 2,..., k} 
corresponds to a vertex, and to each factor fi J we associate an edge between particle J and 
J in the corresponding graph. More explicitly, 


E & 4 (2.6) 


As an example for the k = 3 case, we have the expansion in terms of graphs shown in Figure 
It is clear from the left hand side of (2.5) that there are no multiple edges between two 
vertices, or edges connecting one vertex to itself. 


We also set = - to one. 
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1 2 3 1 3 2 2 1 3 2 
+ e—e—e + e—e—eo + e—e—eo + 1 3 


Figure 1: The graph expansion for the k = 3 case. The second line contains all connected 
graphs. 


This expansion is particularly useful, as the logarithm of the grand canonical partition 
function can be formally given as a sum over connected graphs (see for example Appendix 
A of for a simple derivation): 


o0 k k 
log Ztong = >. T f I] pasion $r) 5 I] fis. (2.7) 
k=1 ` 1 


ge gi" e(I,J) eg 


Here GË! denotes the collection of connected graphs on the set |k] and e(J, J) belongs to the 
graph g if there is an edge between the vertices J and J. Graphs up to four points are shown 
in Figure 


e KT ee 
ra A +e A+al [+6 N+ N 


+ higher-points graphs 


Figure 2: Connected graphs in the Mayer expansion. The coefficients correspond to how many 
different graphs of the given topology belong to Gil, They are obtained from one another by 
relabeling the vertices. 


We stress that (2.7) is an exact relation as a formal power serief] in q. In the limit of 
small € one has 


F 2 
fis = € G(r) + Ole), (2.8) 
3The interesting question of convergence of the Mayer expansion can be addressed in various ways, see 


e.g. (52|. Here, we will not consider this problem or the analog convergence issue for the solution to the 
TBA. 


10 
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where G is defined in (LII). As each factor of f contributes one power of €), the leading 
contribution to the sum in is given only by connected tree graphs TL. Such graphs 
have the minimal number of edges (k — 1) among the connected graphs. The sum of these 
tree graphs can be shown to be convergent. Collecting the powers of €2 we conclude that 


log Zions = 2 A I a6) E M fv+ 04). e 


ge TI! e(I,J)€g 


Graphs oe to five points are shown in Figure[3] Notice that the leading behaviour of the free 
energy (2.9) is at order This fact may be not apes from the definition where for 
each k n leading aontmbuion: is proportional to =. We will see that the similar behaviour 
applies to the more complicated situations ee in the following. 


ee 3 ee ee 
+ 60 e e e © e + 60 e o <L 5 >< + higher-point graphs 


Figure 3: Connected tree graphs up to five-point. The coefficients corresponds to how many 
different graphs of the given topology belong to T”, which are related to each other by rela- 
beling the vertices. 


2.2 Only short range interactions 


Next we study another simplified version of (1.4). Namely, we set D, defined in (2.2), to 
zero and consider 


co k k 2 
Zim = S EL TT aa. (2.10) 


| 2mi € — 
k=0 I=1 A wee N €2 


The main achievement of this subsection is to show that in the limit of small €2 the logarithm 


of (2.10) is given by 


d 
log Zghort = = = = f Buia qQ (o)l + Ole). (2.11) 


We will see that, as opposed to the case of long range interactions considered in the previous 
subsection, the right hand side of (2.11) is not given by summing over tree graphs but by 


11 
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a single local term. As will be explained in details in Section this is the underlying 
mechanism by which the short-range interactions turn a group of particles into a single 
effective particle. The result (2.11) explains why the dilogarithm function appears in the 
TBA. It will also be an essential input in the all order proof presented in Section 
Considering its importance, we will prove in two different ways. 


2.2.1 Dilogarithm from a sum over residues 


In this subsection we will prove (2.11) by direct evaluating (2.10) and (2.11) as a sum over 
residues. We find it instructive to first consider the partition function for the U(1) gauge 
theory, i.e. we evaluate (2.10) for 


1 


Q(x) = G=5-10) @ ee) (2.12) 


compare to (L6). In (2.12) a and € are real and we wrote explicitely the 70 prescription. For 
each k in (2.10), there is only one residue (up to permutation of the integration variables) 
in the upper half plane, compare to the general discussion in Appendix It is given by 
(o1,---,; Pr) = (4,a+€,...,a+(k — 1)eg). It follows that 
“aia, 1 
-ne | ns 2.13 
Short = | ek II & Te ( ) 


One can also directly calculate the one dimensional integral 


| ee) = — = i (aa 


where Q(x) is given in (2.12). Using the identity one concludes the the logarithm of 
(2.13) in the e2 > 0 limit is indeed given by 2.11) with (2.14). 


For the more general U(N) theory, it is easy to classify the poles contributing to (2.10) 
following the same reasoning as in Appendix For fixed k, residues are classified, up to 
permutation of the k particles, by a set of integers {s1, S2,..., SẸ} such that 3, s) = k. 
The corresponding pole is given by 


OP = a+ (h — 1)ez, I= 1,..., S), => Sus (2.15) 


and \ = 1,...,N. As opposed to the full partition function (1.4), for which residues are 
classified by N-tuple of Young tableaux with a total number of boxes equal to k, in the 
simplified integrals (2.10) only Young tableaux with one row, whose length is denoted by 


12 
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sy, contribute. It is straigthforward to calculate and collect all residues, see Appendix [B] to 
obtain 


g? ARA (aaxx + Uy — Dy ee)? 
Se = et = ote 241 
s= D (SR) e 


ian Sl Ol AN D, Ly (aw + 
where 
+( = 1)e2) 1 
Ra TT O ) Os a AU 
s I e + — ) Q (x) POPO (x + €) ) ( ) 


and P(x) = JI wet — a,). The structure of the result (@2.16)-(2.17) represents a simple 
generalization of the U(1) computation (2.13). The e2 — 0 limit of the logarithm of (2.16) 
can be readily obtained using the relation (B.7)). One recognizes that 


Dds T Resio=a [Q (¢0)] + O(€) . (2.18) 


2 St = 


1 
log ZShort = a 


This result coincides with the evaluation of (2.11) by residues. This calculation provides a 
direct proof of (2.11). 


2.2.2 By Mayer expansion and separation of regions 


We now calculate for small €2 by first applying Mayer expansion and then the so-called 
method of expansion by regions. The first step is to decompose the interaction factor as 


Oty 
z= l+ fi, fis := = . 
is — 6 Pig = 6 


2 
£2 (2.19) 


In a similar way as (2.7), the logarithm of (2.10) is then given by 


log Zshort = are TMEL 1) ` I] JIJS (2.20) 


gegt! e(I,J)Eg 


To show that this gives (2.11), we need to evaluate the integrals entering as a 
Laurant series in €2. We are actually interested only in the terms with leading negative 
powers of €). We use a powerful method, usually applied to the evaluation of Feynman 
integrals, called expansion by regions introduced in [49]. The method goes as follows: (1) 
divide the integration domain into regions and expand the integrand in a Taylor series in 
small parameters, (2) extend the integration to the full domain of integration, (3) set to zero 
scaleless integrals. This method turns out to be particularly efficient to calculate the leading 
term for the integrals (2.20). Concerning these contributions, we do not need to apply the 
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Figure 4: The decomposition of the contribution for k = 3 case. Each square represents a 
cluster labeled by corresponding particles in the cluster. 


somewhat subtle point (3) above. A more detailed discussion of this method will be given 
in Appendix 


To identify the set of relevant regions, we note that the Taylor expansion of the interaction 


2 
term F-2 starts at order e2 except for the region in which ¢7, is of order e3. It is then 
Id 


natural to expect that regions are classified as follows. Let Bl"! denotes the set of grouping 
of k labeled particles in clusters, see figure [4] for the k = 3 case. For each grouping we define 
the corresponding region as 


ts ~ €& if I,J are in the same cluster, (2.21) 


lér7| > e if I,J are in different clusters. 


The next step is to Taylor expand the integrand in each region. We denote by T® [...] 
the operation of Taylor expaning ... in the region corresponding to the grouping b € B% 
We have 


(2.22) 


T°) [fr] = O(1) if I, J are in the same cluster, 
‘1’ ~ ) O(e) if I, J are in different clusters, 


and 
Te) Tren] = [J Gy) +... = ra > or (2ra) 


the product on the right hand side runs over the clusters - the grouping b and |Y | denotes 
the number of particles = n cluster Y. From (2.22)-(@.23) it is clear that, for each k, the 
leading contribution to comes from the region in me all ¢; are in the same cluster. 
Indeed, if there were fans or more clusters, the fry factor, with J,J in different clusters, 
would decrease the power of €2, compare to (2.22). An explicit example of this expansion 
of k = 3 is given in appendix [C] As we will shortly see, the leading contribution to (2.20) 
starts at order 1/€, for each k. This fact is not obvious from the form (2.20). 


By further separating the integration over the “center of cluster” coordinate Z as 


1 S 
J Uer (£1, £ y= fa f Tle (3) z) F(@+21,...,+2), (2.24) 
I=1 
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where we applied the change of variables 


on , : 
=z) ts Tr: =T- T, (2.25) 
T= 
one can rewrite as 
Lod 
"naan IŽE jawa oe (2.26) 
2 k= 


where : 7 
1 dé; [1 
- a Sha) D TL tw (2.27) 
I=1 I=1 gc glk e(I,J) Eg 


Notice that in (2.27) all connected graphs, independently on the number of edges, contribute 
to the leading term. The next observation to be made is that 


k k 
1 d 1 
if g is not connected as ft (G X o) I] fi; = 0. (2.28) 
2 T=1 I=1 


(I,J)Eg 


The validity of this statement can be easily argued as follows. For each connected component 
in g we can define its center as the average of the {¢ġr} in that connected component. The 
integrand in does not depend on the distance between the centers. As the integral 
is calculated by residue, it trivially vanishes in this case as, after integrating trivially the 
delta function, the (k — 1)-dimensional residue is absent. The identity implies that in 
(2.27), we can replace the sum over connected graphs with the sum over all graphs. Finally 


using the relation 
DY M f= JL 0+, (2.29) 
geGglk] e(,J)Eg 1<I<J<k 
where GË! is the set of all graphs on [k], we conclude that Jy = Zp, where, using (2.19), 
1 5 dor 1 í o7, 1 kl 
~ F Sai mew aH =a (2.30) 
e71 Ji 2ri k >, euo U = e Ini k2 


Notice that the evaluation of this integral is exact, see for a derivation. Using this result 
we recognize that (2.26) is equal to 


log Zstort = — toe [Bate bo) + (8) = E f SE Liogo] + OF). 231) 


277% 277% 
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2.3 Expansion of the instanton partition function 


We are now ready to consider the full instanton partition function (1.4). We will combine 
the considerations from the previous sections. In Section we learned that the 6 — 0 


limit singles out certain tree level graphs. In Section 2.2] we learned that the factors -— 


produce the effect of combining particles into clusters. To exploit the combination of these 
two mechanisms in the most transparent way we find it convenient to use the so-called 
iterated Mayer expansion advocated at the beginning of our analysis. We start by reviewing 
this expansion. Based on this iterative expansion, by combining the discussion from the two 
previous sections, one obtains a tree graph expansion for the full instanton pationtion as 
given in Section In Appendix [D] we show that the same result can be obtained by 
first applying the method of expansion by regions to the original partition function and then 
exploiting some combinatorics to conclude that only certain connected graphs contribute to 
the free energy. 


2.3.1 Iterated Mayer expansion 


We start with a review of the iterated Mayer expansion [50]. Consider the partition function 


o0 k k 
ee a a ] [ 4¢:A(¢:) I (LEE) 14 Fry = VE Ovi Or), 
k=0 I=1 1<I<J<k 
(2.32) 
where, as in (2.2), we split the pair-wise interaction potential in a short and long range part. 
Next we introduce 


gy = II (1+ ae) l+fry = e Vs" (brs) | (2.33) 
IJEY 
I<J 

1+ Tighe = I] I] e Prat) , a#b (2.34) 


Ia€ Ya EY, 


where Y denotes a set of particles or, in other words, a cluster. Iterated Mayer expansion is 
the statement that the free enery can be expanded as 


oo k k 
y q > 5 f 

log A = E! f J [4% See ce (2.35) 
c= 30 gy, genh” TI 


where 


£ 
Giy, Y} = (u SY. sv.) y I] Ya, Yo (2.36) 


a=1 gE gel e(a,b) Eg 
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-J4 (o) X |] tz, (2.37) 


gEGY elI, Jeg 


The sum in (2.35) is taken over Bi : groupings of k (labeled) particles into @ clusters. As the 
notation may need some time to be digested, in Appendix [B.3] we spell out the definitions 
for k = 2,3. Note that although the cluster here is in a different context, the picture of 
grouping is similar to that used before in expansion by regions, see for example figure [4] for 
k = 3 case. 


2.3.2 Iterated Mayer expansion for the Nekrasov partition function 


We can apply the pial Mayer expansion reviewed in Section to the full partition 
function, compare ) to (4), (2.2). Once this is done we need to evaluate the leading 
contribution for ie € to integrals of the type (2.35). The crucial observation is that the 
Taylor expansion D(¢;;) — 1 ~ €2G(¢,,) is valid everywhere in the domain of pea 4 
This immidiately implies that in the limit of small €2 one can write, see definition (2.36), 


Styr} = & 'II Sy, rD I] I] Gl Omri as (2.38) 


[4] e(a,b) €g la€ Ya 
JETe ( ) wT EY, 


where ... refers to next to leading contributions in €2. We point out that, in analogy with 
(2.9), only connected tree graphs contributes to (2.38). As opposed to (2.9), now they are 
tree graphs on the set [£] of £ clusters rather then the set of k fundamental particles. Notice 
that we did not expand the Sy factor. The integration in (2.35) is still over k variables. 


The next step is to explicitly perform the integration over the distances of particles within 
the same cluster. This turns out to be essentially the same as in Section 2.2.2] For each 
cluster Y, we introduce a “center of cluster” coordinate ¢, as in (2.24). It is clear that the 
fry entering the factors Sy, , see a are independent of the “center of cluster” coordinates 
$a. Using this observation and we write 


[The Sty... Yo} 
I=1 
> JL [| COettrnn)+... (2:39) 


a, duy, I] Q(b, F- Pra) 
gE T}! e(a,b) €g H z y 


=u 


a=1 


Ia EYa 


where 


l- d I 
mom (ee) SM we e 


2 gEGY el,J)eg 


4 Tn the language of the method of expansion by regions we would say that the leading term in the Taylor 
expansion of D(¢;,) — 1 is the same in any region. 
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As extensively discussed in the previous sections, the leading contribution in € of integrals 
of the type (2.39), can be obtained by neglecting the deviation from the center of cluster 
coordinate in the functions Q(x) and G(x). Doing so, (2.39) becomes 


[Thee Ory,,...¥} = + fie. Q™I(G,) Fixa oe I] [Yal [Yo] Gha) + 


ge Te e(a,b) Eg 


(2.41) 
J, was defined in (2.27) and computed in (2.30). For convenience we recall it here 
1 Yal]! 
— d z . 2.42 
Jyul / PYa Oni YAP a) 


In the limit of small €2, the summands in depend only on the sizes n1,...,n¢ > 1 of the 
clusters corresponding to b € B. Converting the sum in from a sum over groupings 
of k particles to a sum over the number of clusters @ and their sizes nı,...,ng produces a 
factor -— aa Assembling the pieces together we finally arrive at the main result 


lim e log Z = 3 >> / J ote) D TI nemla), (248) 


Ni,- Nne>1 g ET}! e(a,b) Eg 


where we replaced ¢, with a. We emphasize once again that the sum has been rearranged 
from a sum over particles, weighted by q, to a sum over clusters. More precisely, for each £ 
there are £ clusters and £ — 1 G-factors. For £ = 1 one immediatly recovers the dilogarithm 


| Erie). (2.44) 


For l = 2 one has 


z | ono) eal) Gls = Qal) = logi- 4Q(6)). (245) 


The result can be visualized diagrammatically in a simple way. For each £ in the 
sum one draws all three graphs on the set {1,2,...,2}. To each node a of the graph is 
associated an integer na > 1. Once the diagrams are drawn, the corresponding integrals can 
be written using the following “Feynann rules” 


e Vertex 2 
nal = £ [Ra on 2.46 
mo] = Sa on.) (2.46) 

e Propagator 
ine 1i = qs G( Gas) - (2.47) 


Here we used a dashed square to indicate that the propagator is associated to the edge 
only. 
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The graphs contributing to @ = 1, 2,3 are given by 


SE a [m] (2.48) 
m Hm (2.49) 


ny 


ny 


L= 2 >. 


n1,n2,n3=1 


ns | (2.50) 


The first two graphs correspond to (2.44) and (2.45) respectively. Notice that each integral 
obtained by applying the Feynman rules, ik = multiplied by an overall ž factor, which 
is explicit in (2.43). Graphs with more vertices, upon dressing the ee with positive 
integers, are as in Figure [B] 


3 Derivation of the TBA 


In this section we show that the logarithm of the instanton partition function in the e2 > 0 
limit, whose structure have been studied in last section, satisfies TBA equations. We will 
first provide a perturbative check to a very high order in the instanton counting parameter 
q, and then present an all order proof based on rewriting the instanton partition function as 
a (0 + 1)-dimensional path integral. 


3.1 Perturbative expansion 


Let us recall the expression of the superpotential coming from the TBA 


w= f [- Foe) (1 - Qe) + Lia aQ(aje*™) ], Ba 
where Fi o 
te) = f SE Gle — y) log (1— aQly)e™) , (3.2) 


compare to (L9) and (@.10). The basic observation is that TBA equation has a natural 
expansion in terms of tree graphs, therefore one can compare it directly with the Mayer 
expansion of the instanton partition function (2.43). 


One writes y and W as 


=X g ole), W= Sd FW. (3.3) 
k=1 k=1 
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y can be solved recursively via the TBA equation, for example up to the first two orders 


alo) = - f Ece- ugo), (34) 
ale) = -3 | Pew-naw?- | EE ce-au- 0E. (85) 


Substituting the expression for y in W and collecting the terms at a given order of q, one 
obtains, here up to order q? 


m =f — (3.6) 
W, = 5 Jm? mi 245 S HH awa - RU) (3.7) 
m = 5 | Eeti | HH eaae- naw) 

+5 / nt Qe) Gle — QW) - AX. (3.8) 


One can see that the expansion of W has a natural interpretation in terms of connected tree 
diagrams, with G(x — y) as propagators and Q(z)’ as vertices. This is the same structure as 
the one obtained starting from the integral representation of the instanton partition function, 
see (2.43). The non-trivial thing to show is that not only the structure is similar but that 
the coefficients actally match. 


By staring at their form, it is not obvious that the Mayer expansion and the TBA 
actually produce the same result. In particular to obtain the contribution of a single graph in 
one needs to take into account a large number of terms coming from the y expansion. 
This procedure can be straightforward implemented on a compute. We have checked up 
to q” order and found perfect agreement with the Mayer-Cluster expansion of the instanton 
partition function. It is remarkable that these two different structures are actually equivalent 
with each other. In the next subsection, we will see the equivalence of the two expansions 
can be proven to all order in q. 


3.2 TBA as saddle-point of a (0 + 1)-dimensional path integral 


As we will shortly review, there is a nice way to rewrite a grand canonical partition function 
of the form (1.4) as a path integral, based on a trick of Polyakov [5I]. Path integrals are 
particularly well suited to be studied in the classical limit A — 0, identified with € in 
this section, where the path integral is obtained by evaluating the classical action on some 
solution to the equations of motion. In the case of the instanton partition function (L.4), 
this limit is more subtle as the two particle potential in the original statistical mechanical 
problem is singular in this limit. Using the considerations from section 2.2] we will be able to 
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disentangle such contributions and derive the TBA in a rather transparent way as a saddle 
point equation of an appropriate action. 


Let us start by reviewing the steps presented in the appendix of [5I]. Consider the 
partition function 


z q" i dor ll V(¢rz) 
Z= a = ers) | (3.9) 
k=0 k! I=1 2ni I<J 
The basic idea is to rewrite the interactions as 
k k 
Lie — Als) [e (3.10) 
I<J I=1 i 
where 5 i 5 5 
A(—) = -5 | tedy Ve (3.11) 
op 2 dp(x) dp(y) 


The operator e^ can be pulled out from the sum over the number of particles k to rewrite 


— AG) dor 6 \| = eA) [Re —o(d0) 
=e aioe [Tee J= Ee 271 


All the information about the interaction in the original gas (8.9) is now encoded in the 
operator A. To translate the construction above into path integral language one observes 
the following equality 


(3.12) 


y0 


AG) Ffy] 


= [De Do] eth actors Hire) 0] F. (3.18) 
y=0 
This can be formally proved by Fourier transformation[Ë]. The measure in [Dy Dp] is defined 
in such a way that if glee = 1 then the resulting path integral should also evaluate to 1. 
Notice ee he relation (8.13) holds also when the operator A is not quadratic. If A is linear 
A= Ly ——., the p eee ae produces a delta function and the resulting path integral 


gives a as i hoali, Indeed in this case e* is just the shift operator. 


From the discussion in Section [2| it is rather clear that, in the study of the instanton 
partition function, we should apply the procedure just described in a sightly modified way. 
In particular we will factor out only the long range part of the potential (2.2), corresponding 
to D, as we will now describe. The first step is to introduce 


ô 1 ô ô ~ 
A(z) a) fa dile =a aaa E O S= (3.15) 
dp 2 dp(x) p(y) 
>The relation (3.13) is a path integral analogue of 
2 dp e xv ip) fẹ dpdy —i i 
AEO, o= f E Ao), = [ G enay), (3.14) 


where f(p) is introduced via Fourier transformation. 
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As above, one uses this operator to rewrite the long range interactions as in (8.10), and 
factor out their contribution from the sum over the instanton number k as follows 


6 ($r) 
Z= AG) Fi] -5i [meem IQ 2k = TI us a _ (3.16) 


=0 e 
P I<J 2 


Notice that we did not expand the short range interaction -— which still appears in 


F|y]. The reader may recognize that we already calculated this re in the €2 — 0 limit 
in Section 2.2} More precisely F[y] coincides with Zgnort defined in (2.10), upon replacing 
Q(x) with Q(x)e~*™. Using the result (2.11) we can thus rewrite 


2771 


Pie] = exp (= f Se Lilaa lo) e] + 016)) (3.17) 


The appearance of the dilogarithm is quite remarkable. The derivation presented makes 
its origin from short range interactions completely transparent. The expression (3.16) with 
(8.17), is well suited for the semiclassical analysis in the € + 0 limit. Using in this 
case, and recalling that Vi (x) = —e,G(x) + O(e) one immediatly obtains 


Z = [Dep ep |= 0.4) +0l@))] (3.18) 


where V[p, p] is defined in and we identified h with e2. In the limit of small e> this 
path integral is calculated by the sadle point methodl4 The superpotential is thus obtained 
from the critical valeus of the functional Y[p,y]. As explained in Section the TBA 
immidiately follows. This complets the derivation. 


3.2.1 Multi-particle interactions 


Let us mention that the above discussion may be generalized to the cases containing “multi- 
particle” interactions, such as 


gh db; qQ(¢1) fis E n EA 
2-5 [Tees oea Tl [I Ontem, B9) 


€ ZJ M=2 h <...<Im 


One can similarly introduce 


An( =) = -1 fer Valers-s20)T] z l (3.20) 


6 Alternatively, one can also integrate out p exactly, as it appears quadratically in the action, before taking 
the e2 > 0 limit. The resulting path integral, in ọ in the limit of small €2, is dominant by the saddle point 
of y, where the saddle point equation is exactly the TBA equation (0.10). 
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In analogy with (8.13), this can be brought to the following form 


Z= [Pope exp - E ae Vne tm) E + J plore) Fly]. 
` aa 


Notice that if m > 2 one can no longer integrate out p exactly. The tri-particle interaction 
is directly related to certain gauge theories. The main example is given by quiver theories 
containing matter in the trifundamental representation of U (2)? [53]. The contour integral 
representation of the instanton partition function for this case can be found in [54]. Notice 
that one has V2 = —e2 G+ O(c), V3 = e2 G3 + O(e). This is the right scaling to produce a 
= factor in the exponent of which is crucial for the semiclassical analysis. 


m=2 


4 Generalization to quiver gauge theories 


Until now the object of study has been the instanton partition function for pure M = 2 
Super-Yang-Mills with gauge group SU(N), see (1.4). This partition function admits a 
natural generalization corresponding to other M = 2 gauge theories. In order to specify 
an M = 2 supersymmetric gauge theory the first step is to fix the following data: a gauge 
group G and a set of irreducible representations R of G corresponding to its matter content. 
Such datum gives rise to a consistent, or UV-complete, theory only if the beta function Øy 
associated to each simple factor of G = [[, Gy satisfies 6, < 0. This requirement puts strong 
constraints on R and Œ. A classification of consistent theories can be found in the recent 
work [55]. In particular, with the exception of G = SU(2)” with R given by trifundamentals 
introduced in [53], each irreducible component of R is charged under at most two Gy. In the 
following we restrict to the case 


G=|[ sum), R= >leon(0,0)+ > 20. (4.1) 


Theories in this family with vanishing beta functions falls into three classes, referred to as 
type I, II, II* in [56]. The matrix c in is then identified with the adjacency matrix of 
the Dynkin diagram of type ADE (type I) and ADE (type II and II*). Any other consistent 
theory in the family can be obtained as appropriate limit of theories of type I, II, II*. 


We fix G and R as in (4.1). We denote respectively by V and E the set of vertices and 
edges of the graph represented by cyw. We denote the number of vertices by V = |V|. The 
instanton partition function depends on the following quantities 


e -deformation parameters €; and €2, we set € = €1 + €2 


e gauge couplings qv, v E€ V 


23 


201708.00365v1 


chinaXiv 


e Coulomb parameters a € Cartan subalgebra of Lie(G), ay, for v € V, A=1,...,M 
e fudametal masses my f, v E€ V, f =1,...,ny 


e bifundamental massed] Me, e E€ E. 


We package masses and Coulomb parameters into polynomials as follows 


Myla) = He = m,f), Rz) = He — yA). (4.2) 
f=1 A=1 


The contour integral representation of the Nekrasov partition function takes the form 


oo k kı kv 
= ESS a qv gauge fund bifund 
Z= > (= -) ky! Eat ky! / II (li Ak (dy) Ak (4v)) I] Že, k (Ose, Pt.) 
“s KV Vv 


ecE 
(4.3) 
where k = J` k, and [dg] = Ky (2ri) dor. As above, V and E are respectively set 
of vertices and edges of the quiver. Moreover, given an edge e, its orientation defines two 
vertices called se (source) and te (target). Using the expressions for z and 2 
given in Appendix [A.3] we rewrite the integrand of as 


I] ti outs) ( I] D(a = a.) ) I] Tl Ul De(se,1 = = (4.4) 


gauge fund 
v » &y 


veV 1<I<J< ky eek LI=1 J=1 
where M,(c) 
x 
v = Py — Ibe Py e ’ 4.5 
Qe) = apes Re- mo) [TT Re +m re (4.5) 
vv! vev! 


x(x + €; + €) 
(a +e) (a+) 

(4.6) 
Note that D is the same as for pure SYM, see (1.5). Moreover, we have the manifest 
symmetry D(x) = D(—x) and D(z) = D,(—2z). 


D(x) = A(2)A(-2), Dele) = [A(me+x) A(m-2))", — Ale) = 


Explanations regarding the origin of this expression are collected in Appendix [A] The 
expression for the contribution to the partition function from bifundamental matter can be 
found in and [58]. For more general gauge groups and representation see e.g. 
and references therein. 


TTo avoid confusion we recall that in our conventions a bifundamental hypermultiplet associated to the 
edge e corresponds to a factor (Os., Dr) + (Oa, Us.) in EI). 
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4.1 TBA for quiver gauge theories 


The partition function in the limit €2 + 0 can be studied in a similar way as for the 
simpler case of pure SU(N) SYM, see (1.4). In particular it can be characterized as solution 
of certain TBA equations given below. The important observation to be made is that as long 
as the bifundametal masses me. are large compared to €2, the pair-wise interaction term in 
the non-ideal gas interpretation of (4.4) splits into two terms with different scales, compare 
to (2.2). The first term e present for each gauge group, gives rise to the dilogarithm 
function in the TBA. The remaining factor in the pair-wise interaction is responsable for the 
kernel term in the TBA non-liear integral equation. 


The resulting TBA is summarized as follows. As in [18], the partition function can be 
written as the critical value of the Yang- Yang functional 


dx d es D 
ES 2ri zA ne 1) GwE) +S far [veelal (wQ,(x)e pv |, 
(4.7) 


where the variations with respect to y,(x) and p(x) give 
p(z) = —log [1 — g Qu(z)e™] , > [te — y)pwly). (4.8) 
The twisted superpotential for the quiver gauge theory is obtained as 
wee — y(x : —py (x 
W = z = (x) log (1 —qmQ (ae ” ( )) + Lig (qv Qulrje ( )) |. (4.9) 
where y(x) satisfy the TBA equation 
dy —ew(y) 
= > zg Cem le — 9) log [1 — qw Qu(yje Po] . (4.10) 


The various functions which contain the data of quiver theory are 


d ete r+ Myw £ — Mw — 41 
Cyl) = dz oe (=) + cw log (==) + cw log ( £ — Mwy )] ’ 
(4.11) 
and Q,(x) defined in (4.5). In (4.11) we defined Myw := Me if the vertices v and w are 
connected by the edge e and mw := 0 otherwise. The expression for the propagator G\w(z) 


given above is obtained as as 


Limites 1 for v=w 
Galt) = Limite s0 E= for v&w (4.12) 
0 otherwise 


Note that in this definition we have Gyw(x) = Gw(—2). 
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5 Conclusion and Discussion 


In this paper we studied the instanton partition functions of four dimensional M = 2 gauge 
theories in a special limit of Q deformation parameters, namely, taking €2 —> 0 but keeping 
€, finite. We show explicitly that the instanton part of the twisted superpotential W, which 
in statistical mechanics language is equal to the free energy, satisfies TBA equations. It is 
thus naturally identified with the Yang-Yang function of some quantum integrable model 
. Based on the proof, we are also able to generalize the correspondence to a general class 
of N = 2 quiver theories. In this case one obtains a set of coupled TBA equations with one 
equation for each node of the quiver. 


The starting point of the derivation is the contour integral representation of the Nekrasov 
partition function. It can be interpreted as the grand canonical partition function of a 
non-ideal gas of particles. Due to the singular structure of the two-particles interaction, 
an interesting effective description of this gas emerges in the €o — 0 limit. In this effective 
description the “new particles” are clusters of the original particles and are subject to certain 
tree-level interaction only. This picture is best obtained by employing the so-called iterated 
Mayer expansion [50]. This structure is nicely captured by the TBA, in particular the 
clustering of particles and their effective interaction is nicely reproduced by the presence 
of the Li, function An all order proof in the instanton counting parameter q of the TBA 
equations is given applying a (0+ 1)-dimensional path integral trick, dating back to Polyakov 
blj. 


In studying the NS limit of the instanton partition function, one has to face the problem 
of evaluating integrals when some parameter, in our case €2, is small. There is a systematic 
way to do so, particularly successful for calculating Feynman integrals, called the method of 
expansion by regions. While our results are obtained also independently of this method, it 
provides a simple way to single out the leading = behavior of the logarithm of the partition 
function. This structure is somewhat reminiscent of the exponentiation of infrared diver- 
gences in gauge theories. It would be interesting to further study the applicability of this 
method to the type of integral considered here. 


In this paper we considered quiver gauge theories involving only U (N) gauge group factors 
and matter in the fundamental/antifundamental or bifundamental representations. It is 
natural to extend this work by deriving TBA equation in the case of the other classical gauge 
groups SO(N), SP(2N) and to other matter content, e.g. symmetric and antisymmetric 
representations of U(N) or trifundamental of U(2)?. Other interesting generalizations are 
the deformation of the partition function by chiral ring operators and the five dimensional 
lift of the partition funtion. The details of these studies will be reported elsewhere. 


One of the most challenging question for the future is to understand how the proposed 
TBA equations emerge from the quantum integrable model point of view. In the case of 


This makes it drastically different from some other TBA obtained by simply applying saddle point 
method [60] [61] [62]. See also [63]. 
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pure SU(N) N = 2 super Yang-Mills, which corresponds to the periodic Toda chain with N 
sites, this has been achieved in [24]. In the derivation a crucial role is played by the so-called 
Baxter Q-operator, which is currently the most powerful and universal tool for determining 
the spectrum of quantum integrable systems. Despite successful applications in many cases, 
see [64| (and also [65]) for recent progress and relevant references, a complete theory of 
Q-operators is still to be established. In particular a systematic construction is currently 
not available in the case for which the representation of the relevant quantum group is not 
of highest-weight type, see and references therein for some progress in this direction. 
Via the NS proposal, quantum integrable systems of this type can be systematically solved 
by gauge theory methods, offering an entirely new perspective on the structure underlying 
integrable models. 


TBA equations of a similar type also appeared in the context of wall crossing effects for 
the BPS spectrum in M = 2 gauge theories [67] [68]. It would be interesting to study the 
connection with our work and find the quantum integrable models corresponding to the TBA 
appearing there. Another connection is to the study of minimal surfaces in AdS [70]. In 
this case, the area of such surfaces is equal, via the AdS/CFT correspondence, to the strong 
coupling limit of the logarithm of null Wilson-loop expectation values in M = 4 SYM. In 
[72], it has been proposed that the exact expectation value can be written as a sum of 
multiple integrals involving some basic building blocks called pentagon transition functions. 
This form is structurally identical to the contour integral representation of the Nekrasov 
partition function. Moreover, the strong coupling limit correspond to the NS limit with 
* ~ €. Techniques of the type presented in this paper can be also used to show, directly 
from this representation of the Wilson loop, that in the strong coupling limit it satisfies the 
TBA equations derived from the analysis of classical strings in AdS [74]. It would be 
very interesting to study the possible connection with our work in this respect. 
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A On the integral representation form of the Nekrasov 
partition function 


In this Appendix we present a brief, non-self-consistent review, together with a collection 
of references, on the origin of the contour integral form of the instanton partition function 


(1.4), B (10). 


The first step is to localize the path integral for the given M = 2 supersymmetric gauge 
theory in the Q-background to configuration of fields satisfying the self-duality equation 
F(A) = 0 and the Dirac equation D4q = 0, see e.g. [75]. Solutions to the self-duality 
equation F',(A) = 0 modulo gauge transformations that are trivial at infinity defines the 
instanton moduli space WM”, The space of solutions to the Dirac equation forms a fiber 
over the instanton moduli space. 


If the gauge group is one of the classical groups U(N), SO(.N), SP(2N) (or a product 
thereof), Mt admits a particular nice description, the ADHM construction [76]. Let us 
review this construction for gauge group U(V) and instanton number k. First introduce the 
linear data 


x = (Bi, Bo,I,J) € Xn, := Hom(V, V) Hom(V, V) Hom(W, V) 6 Hom(V, W), (A.1) 


where V ~ C* and W ~ C™ carry the defining GL(k) and GL(N) action. We refer to (AI) 
as ADHM matrices. On such matrices the group GL(2) x GL(k) x GL(N) acts as follows 


(Ba 1, J) + (go (MEBs) 95", ga 193°, det(M) go J g7’) , (A.2) 
where a, 8 = 1,2, the summation over 8 is understood and 
M e€ GL(2), Ja E GL(N), go E GL(k). (A.3) 
The next step is to introduce the ADHM equations 
ur := IT! — JJ + [B,, Bİ] + [B2, B} =0, (A.4) 


Mc = IJ + [Bi, Bo] =0. (A.5) 


These equations are covariant under the group action given by (A.2). The k-instanton moduli 
space is decribed as 


Myg = {x E€ Xv» such that uele) = 0}/GL(k). (A.6) 


There are two non-trivial steps involved in this definition: (1) impose the ADHM equations, 
(2) divide the resulting space by the action of GL(k). In the following we will see how these 
two steps give rise to the contour integral form of the instanton partition function following 
two different procedures. We anticipate that in both cases the integration over (@1,..., dx) 
in originates from quotienting over GL(k). 
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A.1 D-instanton action and its localization 


There is a natural way to understand the ADHM construction within the language of branes 
[771 [79} [78]. In the simplest example, the starting point is a system of k D(-1) branes and N 
D3-branes . The idea is to look at this system either as the theory leaving on the D3-branes, 
which corresponds to a four dimensional M = 2 gauge theory, or as the d = 0+0 dimensional 
theory on the D(-1) brane. In the latter case one identifies the ADHM equations with D 
and F-flatness conditions of the auxiliary 0-dimensional theory. In this way the instanton 
moduli space Myx is identified with the Higgs branch of the theory on the D-instanton. 


From the point of view of the four dimensional gauge theory the D-instanton partition 
function originates as an integral over the k-instanton moduli space. Following [8], see also 
[75], one can present the integration over the instanton moduli space My., see (A.6), in 
terms of the one over Xy p, see (A.I), and certain “auxiliary fields”. The latter, as emphasized 
after equation (A.6), are introduced in order to implement the ADHM constraints and the 
operation of taking quotient by GL(k). See also [84]. 


The correspoding partition function reads 


Do = 
Zk = | ——~ [DIDY] [DXDH] [DDr] e . A.7 
e= | aay PePIDxDH] [DEn] e (A.7) 
This equation needs some explanation. Measure factors in |...| correspond to BRST doublets 
x Y yY x 
Q X = H 2 Q |H = Oei X i (A.8) 
9 n n 


The variable x corresponds to ADHM matrixes (A1), x = (xr, xc) belong to the same 
space as the ADHM equations (A4), and @ € Hom(V,V). Their BRST partners 
(Y, H,n) have the same transformation properties under GL(N) x GL(k) x GL(2). The 
variable ¢ has no BRST partner, it satisfies Qd = 0. (x, H,¢,¢) are bosonic, (Y, x, n) 
fermionic. The action of tS € Lie(torus) is defined as taking the infinitesimal version of 
the transformation (A.2)) with 


Ga—1 ~ diag(aj,...,an), go—1 ~ diag(¢i,-..-, dx), M-—1 ~ diag(e1, €2). (A.9) 
It is clear from the action of Q that Q? = ts, The action has the schematic form 
So = Q(x-w(x) + tx- H+ tnle, g] + tvor), (A.10) 


where u(x) are the ADHM equations. The parameter t, t’, t” have a fixed value in the original 
problem, but, as the partition function is independent of their values, we can calculate it for 
the most convenient choice. The suffix Q refers to the so-called Q-deformation introduced 
to regularize the infinite volume of moduli space [$] 9]. 
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The variables (7, x, H, $) come into a quartet that can be integrated out in the following 
way [8] £0] [75]. Add to the action a Q-exact term 


6S = sQ(xr- 6) +s Q(xy). (A.11) 


For s > 00, the Hp, yp integration produces the constraints ¢ = 0 = n. Next taking t > oo 
the Hc integration is Gaussian and produces a trivial factor, the yc integration is Gaussian 
(fermionic) and produces the factor 


et, = I] (j-e). (A.12) 


1<I<J<k 


The notation det,,. refers to taking the determinant in the vector space to which yc be- 
longs which is the same as the space of ADHM equations pic, see (A.5). To calculate the 
determinant we also reduce the ¢ integration from G'L(k) to its maximal torus 


dp 
Do > il I] a I] Ory: (A.13) 


I=1 1<I<J<k 


Finally the integration over (x, ~w), is done sending s’ > oo. The (B12, YB») integration 
produces the factor 


1 1 1 
E I oe (A.14) 


and the (I, , Yr), (J, Yz) integration produces the factor 


k N 
1 
im, detr; Q -H II C=: (A.15) 


We have thus SS all integration in (A.7) except for a .--, Ok). Putting together 
(A.12), (A.13), (A.14), (A.15) one recognizes that T given in (A7) coincides with (1.4). This 
ends the discussion of pure a with gauge group U (N). The discussion can be generalized 
to other classical gauge groups and quiver gauge theories, see e.g. [58] and references therein. 


A.2 The instanton partition function from 5d perspective 


As a different picture, we consider the lift of the given NV = 2 four dimensional gauge theory to 
five dimensions. The 2 deformation can be understood as imposing the identification (x, y) ~ 
(ef?x, y+ 3), where y is the coordinate of the circle of compactification and Q the generator 
of an infinitesimal SO(4) rotation. The important fact is that the 5d partition function 
Zé reduces to the twisted Witten index for the supersymmetric quantum mechanics on 
the instanton moduli space My, [12], as fermions have periodic boundary condition (up to 
the twist) along y direction. 
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As descibed in details in [80], by the ADHM description of the instanton moduli space 
(A.6), the 5d index can be ie as 


= i= PADHM e uations) 
gia f d (ei 9) — et (%9) LIU = Pape equations) A.16 
ü W| Vol(T,) -—— (T+) Th I Or i ) [a — WADHM maces) l ) 


The o a equations 20d WADHM matrices are the weights under a a torus action, 
of va is the infinitesimal version, ae = ADHM equations and ADHM ma- 
trices eas The structure of ((A.16) originates from the z step procedure 
E iae after equation (A.6): the term TE — PADHM equations) encodes the constraints 
from the ADHM equation, the integration, and the Vandermonde factor, implements the 
GL(k) quotient. In the following we spell out this formula for pure V = 2 SYM with gauge 
group U(N). 


Consider an element in the maximal torus T C GL(k) x GL(N) x GL(2) given by 

= diag(e!?*,..., ptr), ga = diag(e?®™,..., e884"), M = diag(q,q2) with |q|, |qo| < 
1, compare to (A.2). The eigenvalues of the torus action on the ADHM matrices Ba, I, J 
are 


da ei ors eB (atr) qı Q etb (Praa) , (A.17) 


where a = 1,2, /,J=1,...,k, and à= 1,..., N. One thus obtains 


k 

1 1 
TG oe, P(e 1) P_alq querer) ; A.18 
I] (1 — WADHM matrices) I,J =L P,g(e “AB PT) )P_olu q2 etler) ( ) 

il 1 1 
So kg ee ee =o — n 2 A.19 
lI (1 — WADHM matrices) Bi (1 = q2)* Ag(q2) ( ) 
where 
N 
ao = [[Q—te?™), — Ag(t) = [] Q—te’?*r). (A.20) 
A=1 IZJ 


It is not a coincidence that these equations are in correspondence with (A.15), (A.14). 


Next consider the the contribution from the complex ADHM equations (A.5). Its weight 
under the thorus action gives 


I] (1 — YPADHM equations) = (1 — qı q2)" As(qı q2) : (A.21) 


Compare to (A.12). Finally the Vandermonde is equal to Ag(1). Putting things together 
according to equation (A.16) we obtain 


(=q) Ag(1) Ag(ar Q) 1 
Ple ibti) Plo oe TOP ` A22 
(1 —q1)*(1 — q2) Asla) Asla) 2T] P3(e-*8 tr) Pa(q q2 e-*8 81) ( ) 


The cases of SO(N) and SP(2N) are similar upon modifying ADHM data and related thorus 
action, [80], see also [81]. Matter can also be included in this picture. The limit 6 > 0 can 
be easily taken upon setting qı = e’°*. Notice that the condition |q|, |q2| < 1 translates 
into a small imaginary part for €1,2. 
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A.3 Building blocks for quiver gauge theories 


We now present the recipe to obtain the integrand for the instanton partition function 
corresponding to the class of quiver gauge theories introduced in Section 4] To which we 
refer for the notation used here. From the discussion above it is expected that the integrand 
should be written in terms of the following data: 


(1) an element g= q(€1,2, a, Q, m) € T a SO (4)Lorentz x G sang x Ginst X G'favour, 


(2) a “collection” of weights we, more precisely eigenvalues e’“¢, under this torus action. 


From this data the integrand can be constructed using the following rule 
D mé” I] (we) ™, (A.23) 
2 2 


where ng = +1. 
Let us introduce the following modules for the torus action: 
Wect™, YWec’, Mec”, LeCerR’. (A.24) 


They carry, respectively, the defining representation of GL(N,), GL(ky), GL(ny), SO(4 )Lorenta 
and thus of the maximal torus. We denote the corresponding characters as xwy.m = 
trwv.m(q), for example yy = ye 7, ¢'%. We denote their conjugate modules as Wy, Vy, 

From the description of the tangent space to the instanton moduli space in terms of 
linearized ADHM equations and infinitesimal gauge transformations, see [11], one finds, for 
each gauge group factor 


Ch(TMny x) = xxv — (et — 1)(e'? —l xv xv + e xv xw- (A.25) 


Appling the rule (A.23) to this character one obtains the integrand (1.4). The contribution 
of bifundamental matter is given by 


Ch(O,,O,) = — eb mw Xm XV, — (4 — 1)(e’? — 1) Xy, XV, + eE Xy, X, ] . (A.26) 


Notice the overall minus sign compared to (A.25). Applying the rule (A.23) one obtains 


= ky kw P lows — Myw) Pali + Mw +e) ` 
Ch( lis Ww) = kwNy S ee (A.27) 
vm TT II A(Myw + dy,1 E Qw,J) ʻ ky 


From this expression it follows that the contribution from Ch(C), Oy) + Ch(Ow, Oy) is given 
by 


a E o (A.28) 
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Finally fundamental matter gives the character 


k n 
Ch(O,) = -xvx “SS Am = [] [] r- m). (A.29) 


The expressions above can also be found in [57] 58]. 


A.4 Contour and classification of poles 


The poles contributing to (1.4) correspond to solutions of simple algebraic equations. As we 
will now review following [85], these solutions are classified by N-tuple of Young diagrams 


CI]. Let 


N 
L—€)\("—-€ 
A(t) := | [ (z-a),  E(e):= Ca (A.30) 
z r(x — 6 — €2) 
We want to classify solutions of 
Alzi) []E(ai—2;), for i=1,...,k, (A.31) 
j#t 


modulo permutations of {x;}. Notice that, compare to (1.4), only “half of the poles” of the 
integrand contributes. This is so as a consequence of the choice of contour. It is clear that 
at least one x has to be a zero of A. Without loss of generality we set zı = a), and get k— 1 
equations 
A'(z;) [| Elei- z;), for i=2,...,k, (A.32) 

j>1 

j#t 
where A’(x) := A(x)E(x — ay,). For the same argument as above applied now to (A.32), 
we set A’(x2) = 0. We have to recall that we must keep x; — 2; # 0 and z; — 2; Æ € + & 
otherwise the denominator in E(x) vanishes. Iterating this procedure one finds that solutions 
of modulo permutations are classified by the N-tuple of Young diagrams. For example 
for k = 3 one has 


{a),, Anos ay; } ; {Gig Oras Ads + Ea} , {G45 an, a Ea; an, ot 2€a} ; {a),, an, F €l, ay, F E2} ; 


(A.33) 
where a = 1,2 and the a, are different from each other. 
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B Useful formulas 


B.1 Computing 


In the following we present the explicit calculation of the residue corresponding to the pole 
2.15). First we introduce the notation 


Yı = 4 dis sey Be Ps Yo = {s1,..., 51 + S2}, EE ,Yn = {si t+---+8n_1,...,k}. 
(B.1) 


The next step is to split the interaction term Il zs — Z- in the integrand as 


5 N 5 2 
II 2i z 7 I] I] pi 7 I] I] _ Pity ; (B.2) 


— € — € — € 
1<I<J<k Pry 2 A=1l I<J “IJ 2) \<M Ly Ply 2 


Only the factor in square bracket on the right hand side is singular in the limit in which y 
take the value give by (2.15). For these terms we use the simple result 


SA ob, en 
lim —ay, —(Ih—l1)e = —_, B.3 
Pry +I 7Ox +1 —le2 Il (Prats ( ) 2) Il bry m e SA l ) 
IJEeY\ 


We also use 


a QQ) = 
lim  (@e,41 — ») [I Q(t) = I moa I] cil Naa eee (B.4) 
T=1 


>a € Ie 
Py $1 Xr IEY, 1 + 2 


where 1 
(1) ee ee PO) = — à). B.5 
O Ha = pO) PONe +e)’ (2) He ) eo) 
The result (2.16) immediately follows. 
B.2 Exponentiation of residues 
One has 
love) k oo e f—1 
qr l q 1 d l 0 
ae — 1)e) = aga eL B.6 
d Ha Hy PERE CEE ro PAR A 


as formal series in e€. 
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We present the following generalization of 


“ (q TI, A(ay + (I 1)6)) A Pr (tat (a — Ie? 
> p (Cee) I H UF 


51, SẸN>1 A=1 1<A<X' <N =l Iy=1 
N ow 
1 q 1 d e 0 
as formal series in e. We do not present a proof of this result, but it can be easily checked 
order by order in q. 
B.3 Illustrative Examples for Iterated Mayer expansion 


In the following we spell out the definition (2.36) for k = 2,3. For k = 2 (2.36) reads 


Cio = Sandan = On=—1) De, Guran = fay.) = De- l, (B.8) 
and the sum of this two terms is by the definitions above Fiz. For k = 3, (2.36) reads 


Gra23 = S$f1,2,3} {1,2,3} 5 (B.9) 
Gurp = Su Oa AH > (B.10) 

Gee = faep iae + ian kebeh + Aap fap 
HHRH F238) 1B - (B.11) 


C Expansion by regions for instanton partition func- 
tions 


One recurrent question in this paper is how to evaluate multiple integrals as a Laurent 
expansion in some small parameter, namely e2. There is a powerful method to address 
this question, it is the so called method of expansion by regions, usually applied to the 
evaluation of Feynman integrals [49] [86]. The method can briefly summarized as follows: 
(1) divide the integration domain into regions and expand the integrand in a Taylor series 
in small parameters, (2) extend the integration to the full domain of integration, (3) set the 
scaleless integrals to zero. It has been applied successfully in many situations in the context 
of Feynman integrals, see for a recent study including proof under some assumptions. 
We apply the same idea to the integrals of instanton partition functions. 


It might look surprising that the prescription above works at all. In particular one 
may worry about multiple counting of contributions. Let us describe how this problem is 


°We are actually interested only in the leading + behavior of the logarithm of the partition function. 
y y 8S 8 
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circumvented in a simple example, the k = 2 case in (2.10) 


ra f oe (1) (da) P25. (C.1) 


2nt 20% 12 — & 


The first step is to identify the set of relevant regions. For this, one may note that the 


2 
interaction term zi is different from one only in the region in which ¢ġı2 is of order e2. 


We thus define two separate integration regions as follows 
DS : |o12| ~ O(e2), and pe : loil > €2, (C.2) 


where S and L denote short and long distance respectively, and the whole integration domain 
is given as D = D5 U DŁ. 


The next step is to Taylor expand the integrand in each region. To do so it is convenient 
to perform the change of variables o := elo, $ := $1 — ¢2. In the region Ds one Taylor- 
expands the Q functions for smali[{] ġ, and in region Dr one expands the interaction term 


< 5. Accordingly we introduce the Taylor expansion operators as follows 
Ee 


Ti] = S TUI] = > Qx(b0) o” 1 (C.3) 
TH] = SOT = See») (F). (C4) 

where ; 
j= 5 p3 aoe QEM (g) (C.5) 


After the Taylor expansion, one extends the integration region to be the full region D. 
According to the prescription, one has 


Pa aE, Fe— 5° f TRIN, PL | TN. (C.6) 
k=0 ¥P i=0 YP 


The fact that this equality holds is a very non-trivial statement. In particular one may worry 
about multiple-counting of contributions. One can verify directly upon giving a certain 
prescription for the evaluation of the relevant integrals. Here we provide a formal argument, 
following closely [87], which can be easily generalized to the k-tuple integrals entering the 
instanton partition function. An honest decomposition of will proceed as follows 


r=D fms Df t= ree ( fe mee ft). (C.7) 


10Note that œo will be equal to the poles given by Coulomb parameters a, > €2. 
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In order to show that this is equal to (C.6), one has to show that the terms in the bracket 
give zero. These terms can be written nicely as 


DIOD EDDA OD EDDA fen. cs 


where we use the fact that the double expansion of the integrand is independent of the order 
in which short or long range expansions are applied[4} so that 


TrsT St. (C.9) 


Explicitly, each term in is given as 


do dD ap 
(S,L) al 0 2\k—l 
[ern = 4 S E o E (0.10) 
This is indeed zero since 
fa (¢?)" =0, for integer n. (C.11) 


One can take (C.11) as the analogue of setting “scaleless integrals” to zero in the case of 
Feynman integrals. As in that case, this step is essential in order for the method of expansion 
by regions to work. 


The above considerations seem to be easily generalized to the multi-particle cases. In 
this case one deals with k-folded integrals. It is expected that the double Taylor expansion 
in different regions, compare to (C.9), will produce some factorized term in the form of 
J do (@)". This is set to zero as in (C.II). 


Before proceeding to some concrete examples, let us make some important remarks. 
When applied to the calculation of Feynman integrals, the method of expansion by regions 
presents the characteristic feature of introducing extra divergences in individual integrals. 
These new divergences must be regularized and should cancel after summing over all regions 
in order to reproduce the correct result?2), The situation in the case of the integrals entering 
the instanton partition function is somehow similar. By staring at (C.3)-(C.4) one realizes 
that upon Taylor expansion, one produces arbitrary integer powers of ¢?. These are positive 
powers for TS[7] and negative for T'[/]. The integration along the real line would produce 
higher order UV divergence in the short regions expansion, and higher order IR divergence 
in the long range expansion. As in the case of Feynman integrals these divergences needs to 
be regularized. In the examples below the regularization corresponds to a certain contour 


11 This corresponds to the case of “commuting expansion” in [87]. 

12For example, in the dimensional regularization D = 4 — 2e, this corresponds to higher order + singular- 
ities, corresponding to either UV or IR divergence. The new UV and IR divergences cancel with each other 
from different regions. These divergence may also be understood as from the scaleless integral, where the 
UV and IR divergences canceled within the scaleless integral itself (see examples in ). 
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integral interpretation of the integrals. This prescription works in all examples we considered. 
We recall that the leading term in €, the term we are actually interested in, does not need 
any regularization of this type. 


C.1 Checks of expansion by regions 


In this appendix we check the method of expansion by regions via explicit residue computa- 
tions for two cases. For simplicity, we will consider the U(1) case, namely in the following 
we take 

f(x) = 


phia, 
fie) x—a+e+i10 
where we introduce f(x) for convenience as it does not have any pole inside the contour of 
integration. 


Q(z) = 


=— C.12 
x—a—i0’ ( ) 


Example one 


We first consider the integral (C.I). It is straightforward to compute it by summing over 
the residues of two poles: (¢1, 2) = (a,a+ €2), (a+ €2,a). One has 


Efla) f(a 
F = f(a)f(a+ e2) = Or (C.13) 


k=0 


On the other hand, by the method of expansion by regions, the integral is also given as the 
sum Fs + Fy. One can check that indeed this form reproduces the residue result order by 
order in €2 expansion. It is interesting that Fs contributes only to the terms of odd power 
of €2, while Fy, contributes to the terms of even power. 


Let us spend a few more words on the integrals in Fy, 


f do, doz flor) F(2) 


Qni 2mi (dy — a — i0) (2 — a — 10)(d1 — a) * (C.14) 


The main issue is how to treat the new denominators $7 in the expansion. We will apply the 
following prescription. We first do integration for ¢,. There is a simple pole at ¢; = a +70. 
After taking the residue, the remaining integral for ¢2 has a (2k + 1)th-order pole of ¢2 = a, 
the finally residue is given as 
f(a) f(a) 
2k! 
This prescription works for the next example as well. This provides a strong support for it. 


(C.15) 
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Example two 


We consider an integral which appear at k = 3 order in (2.20) 


r= [TI Law a a (C.16) 
Dai” — 3 ¢?, — E 

As in the example above, one can compute the integral directly by residues. As there is 
no @y, terms in the numerators, also poles for which different @’s are at the same position 
contribute. Let us summarize the result of the residue calculation. The poles of (@1, 2, 3) 
and the corresponding residues are given as follows: 


(a,a,a) : o (C.17) 

(a,a,a + €) + (a,a + €2,a) : a (C.18) 

(a + €,a,a)E4 : r = L (19) 

(a,a + €,a+6) : fla) (at ey" (C.20) 

(a + €,a + €,a) + (at €2,0,a+6€) : La) f(a re lere (C.21) 


One can see that each pole starts at the order 1/e5. However, summing over the residues 
one finds that the first two orders are cancelled, the leading order starts at 1/¢3, which gives 


OPEO SOOO 1) o 
863 3€2 a . 


While for our applications, see , we are interested only in the leading contributions (in 
the normalization used here he is P order =). In the following we will verify that also the 


next to leading contribution is reproduced by the method of expansion by regions. 


Let us compute (C.16) by the method of expansion by regions. According to the discus- 
sion in Section [2.2.2] there are five regions to consider in this case, we have 


F —_— Push + Puth2}.3} + Putt3}.2} + push + PUA ASH . (C23) 


In the first region, there is only one cluster, and one only Taylor expands the Q functions 
around the center of cluster coordinate @, where 


= _ Qı +Q +o 
a 


Q: or := $r — po, (C.24) 


13Note in this case, ¢; = a + €2 should be taken as a double pole, due to ¢12 — €2 and ¢13 — €2 in the 


denominator. 
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and $`} @r = 0. One has 


Sah 3 7 
— do 7 OI = ~ 1 | 
{{1,2,3}} — I I=1 3 2)\ E 
F = [aZ 3) (Q°@) +0x 61+ 067) aaa 
Fai 12 — E2 Pig — €2 
- (C.25) 
The leading term is straightforward to compute. By integrating out ¢;’s one has 
1 f dọ 03 
= fee C.26 
a | O, (C.26) 


which after integrating @o exactly reproduces the leading contribution in (C.22). The re- 
maining terms in the Taylor expansion will contribute to O(e8) order (since ¢7 ~ €3). These 
will constribute to the O(e$) order in (C.22). 


In the next relevant region, there are two clusters. Introduce the center coordinate for 
the two-particle cluster as 


- bite 
= 


Pa : a = {1, 2} , o = 1 — Q2. (C.27) 


and set the coordinate in the one particle cluster 6, = 43, where b = {3}. Doing Taylor 
expansion one has 


FEH — [& doy dd ONAA) + OG, i) (C.28) 


2ri 21 2ri oy (¢? = e2) 


where O(¢, €3) corresponds to higher order terms in the Taylor expansion, which are sub- 
leading in the small €2 expansion. To compute the leading term in this region, one needs to 
evaluate the integral 


dx dy f(x)’ f(y) 


ni Ini (¢ —a—W)*(y —a—W)e—y?’ (C.29) 


where, for simplicity, we considered the case corresponding to the U(1) gauge theory. This 
integral is similar to (C.14). To evaluate it, one can first integrate out x by taking residue 
at the double pole x = a+ i0 then do integration for y. One can also first do integration 
for y and then for x, the result is the same. By summing the leading contributions from 
both regiond4 {{1, 2}, {3}} and {{1,3}, {2}}, one reproduces the next-to-leading order term 
in (C.22). It is clear that the region {{1}, {2}, {3}} where all particles are well separated, 
compare to (C.23), contributes O(e2) to (C22). 


This example provides a non-trivial check for the correctness of the method of expansion 
by regions. In particular, it shows that if one is only interested in the leading contribution, 
the method of expansion by regions turns out to be extremely efficient. 


“The region {{2,3}, {1}} in this case contributes to higher order, since there is no wae factor in the 
23 2 


integral. 
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D Alternative derivation of (2.43) 


In this Appendix we rederive (2.43) by applying the method of expansion by regions before 
taking the logarithm of the e function. This provides a further independent check 
of the result. The validity of the method of expansion by regions is argued in appendix 
here we mainly exploit the combinatorics coming from the labeling of the regions. 


The simple fact on which Mayer expansion is based on is 


J[ < (1+ fr) = `> II i= D II » II fra | > Da) 


1<I<J<k gegil e(I,J)€g be Bll Yeb \geGg¥ell,J)eg 


For convenience we summarize the notation used in the text here: G* (respectively GX), 
where X is a set, denotes the set of graphs (respectively connected graphs) with the elements 
of X as vertices. B* is the set of grouping of the elements of X into clusters. We set 
[A] = {1,2,...,k}. Finally, e(I, J) € g means that there is an edge connecting the vertices 
I and J in the graph g. 


The rewriting together with basic properties of the exponential implies that only 
connected graphs contibutes to the logarithm of the grand canonical partition function. The 
expansion is not particularly convenient when expanding the integrand in a region, 
see Section 2.2.2] There is a refined version of that we will now present. For any fixed 
b € Bll we can rearrange the integrand as 


lu IL @+fy) = ||“ J] G+hy.), (D.2) 


{=l 1<I<J<k Yeb Yı,Y2 €b 
Yi<Yo2 


where 
= I] Ur I] (1+ frs), + fyiy = I] (Lge): (D.3) 
Icey IJeY heYi 
I<J 12€Y2 


The factors Uy, fy,-y, respectively takes into account “interactions” within each cluster Y 
and between pairs of clusters Y1, Y2. Applying the same reasoning as in (D.1) to the right 
hand side of (D.2) we rewrite 


Ho I] G+f)=JT[’ X JJ |S I tw}. @4 


1<I<J<k Yeb be BIN] Y’ Eb’ \geGY' e(¥1,Y2)€9 


The number N, denotes the number of clusters in the grouping b. The sum over Bl) in 
(D.4) is interpreted as a sum over clusters of clusters. 
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Assuming the validity of the method of expansion by regions we have 
Z = a 3 [tr ) [Integrand,] , (D.5) 
be Bik] 


where we define Integrand, as the r-h.s. of (D.4). The symbol T® [...] means the following: 
Taylor expand its argument in the region labeled by b. Based purely on combinatorics one 
rewrites 


log Z = =e 2 ees Il’ dS |i ea (D.6) 


be Blk Yeb geg! e(¥1,¥2)€ 


The explicit Taylor expansion gives 


b) [iy] = (2%) Ti 


oy := Ea P $r, (D7 


I, JEY p -gt IceY 
Tred 
T®) ben el = č |Y;| Y>] G(oy, E y) gers. (D.8) 


From one concludes immediatly that only connected tree graphs contributes to the sum 
in (D6). This produces a factor of €}°~'. Next, for each cluster we separate the integration 
variables into center of the cluster and distances between particles in the same cluster. The 
latter contribution facorizes from to produce 


dor, ys 1 
F af I a (a La), Uoo oap 
I<J 


We have thus reproduced the result (2.41). 
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